so (4) symmetry in the relativistic hydrogen atom 
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We show that the relativistic hydrogen atom possesses an SO (4) symmetry by introducing a 
kind of pseudo-spin vector operator. The same S0(4) symmetry is still preserved in the relativistic 
quantum system in presence of an U(l) monopolar vector potential as well as a nonabelian vector 
potential. Lamb shift and SO(4) symmetry breaking are also discussed. 
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Symmetry principle is one of the cornerstones of mod- 
ern physics. It has been playing a more and more sig- 
nificant role in theoretical physics since the early twen- 
tieth century, when Einstein first put it as the primary 
feature of nature that constrains the allowable dynami- 
cal laws The Einstein's profound change of attitude 
on symmetry principle has made a great progress in the 
study of symmetry. In the latter half of the twentieth 
century it has become the most dominant concept in the 
exploration and formulation of the fundamental laws of 
physics, and nowadays it serves as a guiding principle in 
the search for further unification theory. Dynamical sym- 
metries are prevalent in many important physical models. 
For instance, in the nonrelativistic quantum mechanics, a 
three-dimensional harmonic oscillator has an U(3) sym- 
metry, a three-dimensional hydrogen atom has an SO (4) 
symmetry and the Haldane-Shastry model, which de- 
scribes the one-dimensional long-range spin-interaction 
chain, has a Yangian symmetry [3|. 

As is well-known that the nonrelativistic quantum me- 
chanics is an approximate theory of the relativistic one. 
This gives rise to a fundamental open question: Sup- 
pose a nonrelativistic quantum system possesses a certain 
dynamical symmetry, when its corresponding relativistic 
quantum mechanical version is taken into account, will 
the same symmetry still reside in the system? Harmonic 
oscillator and hydrogen atom are the two simplest pro- 
totype models in quantum physics. Dirac himself has 
introduced a kind of relativistic version for the quan- 
tum mechanical harmonic oscillator with the Hamilto- 
nian H = a ■ {p — i[3Mujr) + 13 M, which now known as 
the Dirac oscillator 0. However, so far the full symme- 
try of the Dirac oscillator has not yet been clear. Very 
recently, Ginocchio has made a remarkable progress [4'| 
by showing that U(3) symmetry does reside in a kind of 
relativistic harmonic oscillator, whose Dirac Hamiltonian 
reads H = d-p + f3M + {1 + (3)Muj'^r'^/2, where a, /3 are 
the Dirac matrices, pis the three-dimensional linear mo- 
mentum, r is the spatial coordinate, r its magnitude, M 
is the mass, lo is the frequency, and the velocity of light 
and the Planck constant have been set equal to unity, 
c — h ^ 1. This fact also suggests that Ginocchio's ver- 
sion of the relativistic quantum harmonic oscillator may 



be a more natural extension than the Dirac oscillator. 

Whether the relativistic hydrogen atom (RHA) has an 
SO (4) symmetry is still open? The purpose of this Brief 
Report is threefold. First, we show that there is indeed 
an SO (4) symmetry in the usual relativistic hydrogen 
atom, by introducing a kind of pseudo-spin vector op- 
erator. Second, we illustrate that the same SO (4) sym- 
metry is still preserved in the relativistic quantum sys- 
tem in presence of an U(l) monopolar vector potential. 
Third, we find that the relativistic hydrogen atom still 
possesses an SO (4) symmetry if some kinds of appro- 
priate nonabelian vector potentials are presented. This 
reflects that the hydrogen atom (or the Kepler system) 
is a highly symmetric system, whatever in the levels of 
classical mechanics, quantum mechanics or even the rel- 
ativistic quantum theory. 

80(4) symmetry in the usual RHA. The Dirac Hamil- 
tonian for a relativistic hydrogen atom reads 



Hrha = d-p + PM , 
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a is the vector of 
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Pauli matrices, 1 is the 2x2 identity matrix, a — e^ 
the fine structure constant, and e the electric charge. Its 
energy spectra is given by the Sommerfeld formula Q 
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where n is the radial quantum number, k = ±(j -I- 1/2) 
are eigenvalues of the Dirac's operator K — /3(I] • L + 1) 
with _ftr^ = + 1/4, J = L + 5 is the total angular mo- 
mentum, L — r y. p the orbital angular momentum, and 
S = S/2 the spin-1/2 angular momentum. In the nonrel- 
ativistic limit, the Hamiltonian ([T]) reduces to the usual 
nonrelativistic hydrogen atom Hnrha = p^ /2M — a/r, 
and the Sommerfeld formula reduces correspondingly to 
the Bohr formula = -Ma'^/2n'^ k. E- Mc^. 

The Hamiltonian Hnrha commutes with L and the 
well-known Pauli- Runge-Lentz vector which form 
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a dynamical symmetry group of S0(4). Evidently, the 
Bohr's energy formula depends only on the principal 
quantum number n, and it has n^-fold degeneracies due 
to the S0(4) symmetry. However, the Sommerfeld's en- 
ergy formula is j-dependent, for fixed n and j, the en- 
ergy has only 2(2j -I- l)-fold degeneracies for n 7^ |k| and 
(2j -I- l)-fold degeneracies for n = |«;|, this does not sup- 
port apparently that the RHA still possesses an SO (4) 
symmetry, since any reduction or elimination of degener- 
acy usually implies the broken of the symmetry. Never- 
theless, the Hilbert space of H^ha is larger than that of 
Hnrha by Considering the additional intrinsic spin space. 
Therefore it is still possible to restore an SO (4) symme- 
try for Hrha through combining properly the operators 
of r, p and E. 

If this is the case, the question is, what are the six rela- 
tivistic generators? What we need to do first is to find out 
six linear independent operators (i.e., six integrals of mo- 
tion for the RHA) that all commute with H^ha , and then 
arrange them to be six generators of SO (4) group. Up to 
now, people have known five of them. The first three are 
three components of the total angular momentum oper- 
ator J. The fourth is the Dirac's operator K mentioned 
above. The fifth integral of motion was discovered by 
Johnson and Lippmann Q in the year of 1950, which 
now called the Johnson-Lippmann (JL) operator. Such 
a famous discovery has stirred a great furor, and many 
people have been attracted in this problem 0, [l^, 11 1. 
The JL operator reads 
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with its square is 
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7^E, and f = f/r is a unit 



vector. The physical significance of the JL operator in 
the nonrelativistic limit is nothing but a projection of the 
Pauli-Runge-Lentz vector on the spin angular momentum 
vector [1, |Tl[. 

The commutation relations among Hrha and these 
five conserved quantities are [i?r/!.a,t^ = [HrhaiK] — 
[Hrha,D] ^ 0, [J,K] = [J,D] = 0, and remarkably 
{K, D} = KD + DK = 0, namely K and D are anti- 
commutative. As usual, the simultaneous eigenfunctions 
of {Hrha.J^.J^} are twofold Krammer's degeneracies, 
i.e., 
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with Hrha\i'tjm,) = ^^l^rym,)- J^W' 

-J to j. Here M — Jq°° dr[f^{r) 
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„jm,/ ■■-ji-rnjm.n and TUj ruus from 
Here Af = Jo°^ dr[p{r) + g^ir)] is the 
normalized coefficient, /(r) and g{r) are real func- 
V l + m + lYi^{e,ip) 
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spherical harmonics, and {a 
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jm, vjm,- The eigenstates \i^%^^{r)) are dis- 



tinguished by the Dirac's operator as i^lV'njmj (^)) = 
±|k||i/)^^„j . (r)). The existence of the JL operator is the 
direct reason that causing the twofold Krammer's degen- 
Later on we shall show this fact from the 



eracies 

viewpoint of the pseudo-spin operators. 

The anti-commutativity between operators K and D 
motivates us to introduce the pseudo-spin vector operator 
T= (ri,T2,T3) = (ri,T2,r3)/2, where 
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The operators T3 and t\ are defined by rescaling K and 
D such that ~ t\ The operator T2 is defined by 
the commutator [ra, ri] = 2iT2, or T2 = itit^. The vector 
operator r plays a similar role as the Pauli matrices vec- 
tor a. If \ipnjmj (^)) is eigenstate of Hrha and T3, then 
'^il^njm (^)) is ^iso an eigenstate of Hrha and T3 because 
of T1T3T1 — — T3. This is the reason causing the twofold 
Krammer's degeneracies. 

It is easy to show that [Ti,Tj] = ieijkTk and = 
i(i + 1). The vector operator T has a property like 
spin-1/2, yet it is not a spin, because it contains r and 
p, consequently we call it a pseudo-spin- 1/2 vector oper- 
ator. One may have [Ji,Tj] = 0, in other words, J and 
T are two independent angular momentum vectors that 
commute with the Hamiltonian Hhra- Therefore, after 
making the following simple linear combinations 



I = J + T, R = J~T, 



(7) 



one arrives at an S0(4) algebraic relation: = 
i^ijklk, [h,Rj] — i^ijkRk, [Ri,Rj] = ieijklk- Since 
[Hrha,I\ = [Hrha,R] = 0, this cuds the finding of S0(4) 
dynamical symmetry in the usual RHA. 

Full energy spectra have been derived by the well- 
known ladder-operator procedure The symmetry 
involved in the system is helpful for us to obtain the 
energy spectra. For example, in Eq. since 
is positively defined, its minimal eigenvalue is zero, 
then one can get precisely the ground state energy 

' — 1/2 

of hydrogen atom [11,] as i?|„=K=i — M (l — a^) 
Furthermore, with the aid of the pseudo-spin-1/2 
vector operator, one can establish an elegant formula 
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between energy spectra E and the integrals of func- 
tions /(r) and g{r). More explicitly, let us denote 
|±) = I'fpnjmX which are eigenstates of T3 with 
eigenvalues equal to ±1/2. r+ = Ti + iT2 is the rais- 
ing operator T+|-) = |+), hence = 1. 
It is easy to check that T+|— ) = "^^fK) = 

-^(i?-M)(j + i)/(r)0f^ 
Thus, (+lTj-) = ^(l-^/r 
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1. Let dr[2f{r)g{r)]/Af = b, we then obtain an 

equation 1 -I- VI?^ = {j + l/2)6/a. By using Eq. ([4]), we 
obtain the desired formula 
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The exact solutions of f{r) and g{r) are [l3j: 



/(r) VM + .B[-nF(l -~n,2v+ 1, p) -|- (AfaA + k) x 

F{-n,2i^+l,p)] p^-^e-P^^, 
g{r) = \/M - E[-r~iF{l - n,2v + 1, p) - {MaX + k) x 

F{~fi,2u+l,p)]p^-'e-p/\ (9) 



the monopolar Sommerfeld formula 
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where k' = ±\/ (j' + 1/2 -|- q){j' + 1/2 — q) are eigen- 
values of the Dirac's operator K', and n' — -I- 
1 + 2 + |k'|, • • • . When q = 0, all the extended 
monopole-dependent operators and relations reduce to 
the usual ones in RHA. Furthermore, in the nonrelativis- 
tic limit, the Hamiltonian (fTU)) reduces to the nonrel- 
ativistic monopolar- hydrogen atom H'^^f^^ = -k'^ /2M + 
q^ /2Mr'^ — a/r. Such a quantum system still has an 
SO (4) symmetry due to the monopole-dependent orbital 
angular momentum vector and the monopole-dependent 
Pauli-Runge-Lentz vector [17 1. 

SO (4) symmetry in the RHA with a nonabelian vector 
potential. Let us add a nonabelian vector potential A — 
iW{r)H X r to RHA, where W{r) ia an arbitrary real 
function of r. Then the Dirac Hamiltonian reads 



Hrha ^a-{p- eA)+(3m- 



(12) 



where n — n — \k\, ly = \^K^ — a^, p — 2r/MaX, A = 
l/\/M2 - £;2. Substituting Eq. ^ into Eq. one 
can verify directly the correctness of Eq. fH]). 

S0(4) symmetry in the RHA with an U(l) monopole. 
In order to find the reason for the existence of the small- 
est electric charge, Dirac published a paper in 1931 which 
started the subject of magnetic monopoles [IJI- After 
1931, the theory of monopoles has been studied exten- 
sively in many literatures, such as [isl [l^ . The Dirac 
Hamiltonian of the U(l)-monopolar RHA reads 



Kha = a-TT + mf3 , 

r 



(10) 



where if = p — eA, A is the vector potential of an 
U(l) Wu-Yang monopole with strength g [l^, and B — 
X A — g^ \s, the magnetic field satisfying the Coulomb 
gauge \I ■ A — 0. Similarly, one finds that H'^^^^ pos- 
sesses an SO (4) symmetry, the corresponding generators 
axe /' = J' -h f ', B! = J' - f ', where J' = V + S, 
L' = fx fr—qf is the monopole-dependent orbital angular 
momentum vector, q — eg — ^ x integer is the magnetic 

charge, f' = {D'/^/d^, iD' K' D'^ K''\ K' /^fx^) /2, 
K' = /3[S-(fx7r)-hl] = l3[t-L'+qt-f+l] is the monopole- 

dependent Dirac's operator, K'^ — J' + 1/A — q^, 
D' = 7^0 -f - {iK' / Ma)j^ {H'^^^ - M (3) is the monopole- 
dependent JL operator, and D'^ = 1-1- {H'^f^^/ AP — 
i){K' /a^). It is worthy to mention that D' and D' 
keep the same structures of D and D^. 

The SO (4) symmetry of H'rha can be checked di- 
rectly. Correspondingly, its energy spectra is given by 



Similarly, H^f^^ possesses an SO (4) symmetry with the 
corresponding generators are / = J + T , R = J — T . 
Here J is the usual total angular momentum oper- 

{d" 



ator, T 



'D"^,iD" K/VD""K^,K/VK^)/2, 
D" ^-f'^a-f- {iK/Ma)-f^{H'^^^ ~ MP) is the JL oper- 
ator in presence of the nonabelian vector potential, and 
d"^ ^ 1 -f- {H"lf^^/AP - l){K'^/a'^). Also the operators 

D and D share the same forms of D and D . The 
nonabelian vector potential satisfies the Coulomb gauge 
V • 0, from = (l/2)e,jfc(9jA - dkAj + [A, A]) 
one has the "magnetic" field as B = i[2W{r)~rW' (r)]T,+ 
i[2W'^{r) + W' {r)/r]{T, ■ f)f. Interestingly, if we choose 
W{r) = -l/r2, then the "magnetic" fields = (-4i/r2)I] 
is proportional to the spin vector. 

Lamb shift and SO (4) symmetry breaking. The spec- 
trum formula ^ indicates that the levels having the 
same n and j values should be degenerate, for instance, 
the 2Si/2 and 2Pi/2 levels should share the same energy. 
However, In 1947 Lamb and Retherford made an elabo- 
rate experiment and it showed that the 2Pi/2 energy level 
was depressed about 1.057 x 10^ Hz below the 2Si/2 en- 
ergy level This effect is now called the Lamb shift. 
It cannot be explained in the framework of the ordinary 
quantum mechanics and gives rise to the birth of the 
quantum electrodynamics (QED). In fact, the success- 
ful calculation of these small quantum corrections to the 
Dirac energy levels was one of the remarkable achieve- 
ments of quantum field theory. From the viewpoint of 
QED, Lamb shift is caused by the vacuum polarization 
and vertex corrections |l9|, and the effective potential 
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FIG. 1: (Color online) Low-lying energy levels of Hrha and 
H'^f^^ with q = 1/2 (not drawn to scale). The green lines 
denote the energy spectra E in Eq. ©, while the red lines 
denote the energy spectra E' in Eq. . The energy levels of 
HL 

ha are depressed below those of Hrha, and such depressions 
fade away as n, j tend to infinite. For n = 2, j = 1/2, the 
depression caused by monopole (« 10^'' Hz) is much bigger 
than that of Lamb shift (?a 10^ Hz). 
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the second term represents the interaction between spin 
and orbital angular momentum. Obviously, the corrected 
Hamiltonian Hrha = Hrha+'^VLamb commutes only with 
J and K, thus the SO (4) symmetry is broken. 

In Fig.l, we have plotted some low- lying energy levels 
of Hrha and H^f^^ with q = 1/2 (not drawn to scale). 
The green lines denote the energy spectra E in Eq. 
while the red lines denote the energy spectra E' in Eq. 
(fTTj) . One finds that E'^^ij, are depressed below Enj. For 
instance, let A^^- = Enj — E'^,j, denote the depressions, 
then for q — 1/2 one has the depression for the ground 

state as ^\^i^2 ^ 1-098 x 10^^ Hz, and that of the first 

excited state as ^\^i/2 — 1-225 x lO^"* Hz. Actually, such 
depressions fade away as n, j become larger, e.g., the 
depression of the levels 3D^/2 (1.294 x 10^^ Hz) is only 
about one thousandth of that of 15i/2 (1.098 x 10^^ Hz). 
For fixed n, j, the depressions caused by monopole are 



much bigger than those of Lamb shift. 

In conclusion, we have shown that the relativistic hy- 
drogen atom possesses an SO (4) symmetry by introduc- 
ing a kind of pseudo-spin vector operator. The same 
SO (4) symmetry is still preserved in the relativistic quan- 
tum system in presence of an U(l) monopolar vector po- 
tential as well as a nonabelian vector potential. When 
the effect of Lamb shift is taken into account, the S0(4) 
symmetry in the quantum system is broken. 
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